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ABSTRACT 

Within the framework of the operator product expansion (OPE) and the renor- 
malization group equation (RGE), we show that the temperature and chemical 
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'sj" ■ potential dependence of the zeroth moment of a spectral function (SF) in an 

• asymptotically free theory is completely determined by the one-loop structure 

. of the theory. This exact result constrains the qualitative shape of SF's, and 

I implies striking phenomenological effects near phase transitions. 

I Our present understanding of QCD at finite temperature (T) and baryon density 

(or chemical potential /x) is mainly limited in the Euclidean realm, due to the lack of 
Oh. non-perturbative and systematic calculating tools directly in the Minkowski space. The 

r-| ! typical methods, with QCD Lagrangian as the starting point, are the OPE and lattice 

simulation. Because of these two formulations are intrinsically Euclidean, only static 
. quantities are conveniently studied. In order to gain dynamical informations, which are 

rS I more accessible experimentally, the analytic structure implemented through dispersion 

' relations often have to be invoked within the theory of linear response. 

The real-time linear response to an external source coupled to a renormalized 
current J{x) is given by the retarded correlator: 

K{x-T,fi) = 9{xo){[J{x),Jm)T,,, (1) 

where the average is on the grand canonical ensemble specified by {T,fj,). Disregarding 
possible subtraction terms, the following dispersion relation for the frequency depen- 
dence of the retarded correlator can be written: 

K{uj,k;T,fi) = du — — - — — . (2) 
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For convenience, we discuss only the uniform limit (k = 0). Upon analytic continua- 
tion, LV ^ iQ, the dispersion relation in principle connects the correlator in Euclidean 
region to the SF, which embodies all the real-time information. However, given only 
approximate knowledge of K{iQ; T, /j), either via the OPE or lattice calculations, it is 
extremely difficult technically to "invert" the dispersion relation. Only when we have 
enough understanding of the SF involved this "inverting" can be made practical. In 
fact, the success of the so-called QCD sum rule approach and analyzing lattice data 
under many situations rests, to some extend, fortuitously on the successful parameter- 
ization of SF's. 

On the other hand, the lack of adequate understanding of SF's at finite (T, /i) 
severely limits our chances to extract physical information using these approaches. It is 
likely that physical results at finite (T. //,) can be strongly biased, if one naively assumes 
the same parameterizations of the SF's that work at zero {T,fi). The purpose of this 
work is to derive exact sum rules that constrain the variation of SF's with (T, /x). The 
derivation, based on the OPE and the RGE, has close analogies to the analysis of deep 
inelastic lepton scattering experiments. In addition, we apply these sum rules to the 
chiral phase transition, and demonstrate that SF's in some channels are drastically 
modified compared both to their zero (T, /i) and perturbative shape. 

In an asymptotically free theory, the OPE yields, e.g. in the MS scheme, the 
large- asymptotic expansion 

K{tQ; T, fi, n) ~ Ko{iQ, + ^ C^iQ', q^k), «:) ([0„],)t,^ , (3) 

n 

where g^{K), [OrJ^'s and C„'s are, respectively, the coupling constant, the renormalized 
composite operators and their corresponding Wilson coefficients at the subtraction 
mass scale k. It is important to notice that the information of the ensemble average is 
encoded in the matrix elements of the composite operators, while the Wilson coefficients 
and Kq are independent of T and /i. Although the matrix elements ([0„]k;)t,^ cannot 
be determined perturbativcly, the Q^-dependence of Wilson coefficients C„ is dictated 
by the renormalization group equation and given by 



Cn(Q^5'(/^),/^) 



g'(Q) 



g'(^) 



(27J-Tn)/2f. 



{l + 0{g'm}, (4) 



where dn is the canonical dimension of the operator On minus the dimension of K and 
Cn{g'^{Q)) is calculable perturbatively. The pure numbers 7j {i — J, n) and b are related 
to the anomalous dimensions oi J, On and to the /3-f unction as follows 

= -7i/ + Oig') , and /3 = -bg^ + 0{g') . (5) 

To study the dependence of K on (T, ^) we only need to consider the difference 
Ak(iQ) = k(iQ;T,ii) - k(iQ;T',iJ.') and 



«)^r*'-^. (6) 
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where Ap{u) = p{u; T, fj,) —p{u; T', p'). This subtraction is crucial to remove Ko^iQ, k), 
which contains all the terms not suppressed by a power of and also to make 

AK{iQ) independent of the renormalization point k. Finite masses give corrections of 
order m'^{Q)/Q'^, with m?{Q) that runs logarithmically and hence can be ignored, if 
we are only interested in the lowest moment of the subtracted SF. 

At this point we have expressed the left-hand side of Eq. as an asymptotic 
expansion of the form: 

n,v=0 ^ 

where A([0„]k) denotes the difference between the expectation values of [Onlre in the 
ensembles specified by (T, /i) and (T', p') and the exponent r]n and the Q^-independent 
coefficients c^n\^) known perturbatively. 

We proceed by making an analogous asymptotic expansion of /S.p{u): 

M-)-Ei$ii5^f:^H/Hr- (8) 

For notational clarity, we have ignored exponentially suppressed terms and the fact 
that there can be more than one r]^ and ^„ for each n. We then obtain the sum rules 
by imposing that the asymptotic expansion of the right-hand side of Eq. (|^), which we 
get by inserting Eq. (|D in the^dispersion integral, matches the left-hand side obtained 
by the OPE. In a long paperEJ we shall present those technical details that make the 
procedure sketched above rigorous. At the moment we only need to consider the leading 
terms in the expansions of, respectively, the left-hand side and the right-hand side of 
Eq. (I): 

[9'm'"^^^[cf{^) + c^:\^)9'm , (9) 



Q 



A-p+^^[g^m'-' 



(10) 



In Eq. ([l^), if ^„ > 1, Ap can be shown to be equal to the zeroth moment of the 
subtracted SF; note that the zeroth moment of a function whose asymptotic expansion 
is Eq. (H) is infinite, if ^„ < 1, since J^dx {log x)^^" = oo. 

For the sake of concreteness, let us examine the consequences of matching Eq. 
and Eq. (p!oD in the case we are concerned with, i.e. n = 1 and di = 2. If the OPE 
calculation produces = 0, then ^„ must be an integer greater than one, the zeroth 
moment exists and is given by 

oo , , 

du' Apiu) = cf\K)A{[0,l) . (11) 







If ?7„ > 0, then = ^ + Vn, the zeroth moment is again finite and equal to zero: 

/ du'^Ap{u) = 0. (12) 
Jo 
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Our main results, Eq. ([TT]) and Eq. ([12|) , can be expressed in physical terms as 
follows. The zeroth moment of a SF for a current J whose OPE expansion yields r/n > 
is independent of T and /x, while the same moment for a current with rj^ = Q changes 
with T and /i proportionally to the corresponding change (s) of the condensate (s) of 
leading operator (s). When ^„ < 1 or ^„ = 1 in Eq. (||) terms such as [g'^{Q)]~^ or 
ln[5f^(Q)] in Eq. (|10D would be produced. Hence the appearance of the inverse power 
of the g'^iQ) or ln[(7^((5)] in the OPE series is a indication that the zeroth moment of 
the subtracted SF is infinite. 

At this point several more general comments are appropriate: 1) It is essential to 
take into account the QCD logarithmic corrections, for the logarithmic corrections not 
only dictate whether Ap satisfy Eq.(|ri|) or Eq . (|T^) , but also control the very existence 
of Ap. 2) The derivation of sum rules for higher moments of the SF requires the 
complete cancelation of all the lower dimensional operator terms, not just the leading 
g^{Q) terms; in particular, we also need current quark mass corrections to the Wilson 
coefficients. An appropriate subtraction is the prerequisite for the convergence of higher 
moments. 3) The (T, p)-dependent part of the leading condensate appearing in Eq. { ^A\j 
does not suffer from the infrared renormalon ambiguity, because only the perturbative 
term Kq can generate contributions to the leading condensate that are dependent on 
the prescription used to regularize these renormalons. But Kq is independent of T and p 
and any prescription dependence cancels out when we make the subtraction in Eq. (j^). 
On the contrary, unless we generalize Eq. (j^) and make other subtractions, sum rules 
that involve non-leading condensates are, in principle, ambiguous. 4) We have explicitly 
verified the correctness of our results in a soluble model,EJ the Gross-Neveu model in 
the large-iV limit, where we obtain exactly all the relevant quantities, such as SF's 
at arbitrary (T, p), Wilson coefficients, /?- and F-functions in pseudoscalar and vector 
channels. 

Now let us specialize to QCD and consider four correlation functions: two in- 
volving the non-conserved scalar Jg = ipip and pseudoscalar Jp = ip'j^i' currents 
ilJs = iJp =_l/47r^) and two involving the conserved vector Jy = ip'j^i'ip and axial- 
vector J A = currents {■jy = 7a = 0). Since the leading operators (dimension 
four) have non-positive anomalous dimensions, the two non-conserved currents have 
Vn > (27j — '~fn)/'^b > and Eq. (p^) applies, i.e. the zeroth moments of their SF's are 
independent of T and p. On the other hand, the two conserved currents have rin = 
and a generalization of Eq. (pA]) applies.El Since there are three dimension-four operators 
with zero anomalous dimension, the sum rules for the vector and axial-vector currents 
are 

/ du'Apiu) = aA{[mi,^]) + + 8A([goo]) , (13) 

JO ivr 

where 6fj_i, is the traceless stress tensor and a = 6 for vector and a = — 10 for axial 
channels respectively. This exact sum rule should not be contaminated explicitly by in- 
stantons, although the value of the condensates certainly have instanton contributions. 
The reason is that the instanton singularities in the Borel-plane are located on the 
positive axis starting at Svr^, and, therefore, contribute to correlation functions only 
with higher order terms of 1/Q^ in the OPE series. 

Finally, let us discuss some of the phenomenological consequences of these exact 
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sum rules. In the pseudoscalar channel, Ap = implies that, in the broken- chiral- 
symmetry phase, the change of the pion pole induced by T or /i is exactly compensated 
by a corresponding change of the continuum part of the SF. Next let us consider the 
scalar correlation function at = 0, the chiral susceptibility, 

x(T,/i) = Jd'x9{xo){[Js{x),Jsm)T,, = l^du' ^^""'J'^K (14) 

which diverges when (T, /i) approaches the phase boundary, provided the chiral restora- 
tion is a continuous transition. The divergence of the chiral susceptibility near phase 
transition only can be produced in Eq.(|l^ by singularities very close to the origin, 
when the exact sum rule Ap = is simultaneously taken into account. Thus the spec- 
tral function would have to have a vanishing threshold (since there is no massless poles 
in the chirally symmetric phase) and develop a strong peak right above the threshold, 
when (T, p) is very close to the phase boundary. Because in the chirally symmetric 
phase the pseudoscalar and scalar channels are degenerate, the same would also hap- 
pen to the pseudoscalar SF. This strong peak in the pseudoscalar and scalar SF's, 
which is intimately connected with the critical phenomenon of diverging susceptibility 
and the correlation length near the phase transition, can be interpreted as some kind 
of quasi-particle, thus confirming the qualitative picture, originally proposed in the 
context of the Nambu-Jona-Lasinio model,i of the appearance of soft modes near the 
chiral phase transition. Similar situation, though less drastic, can be argued to occur in 
vector and axial channels, based on the facts that A ([6*00]), A{[miljip]) and AdagG"^]) 
behave smoothly across the critical line and the rapid increase of the so-called baryon 
number susceptibility at the chiral restoration point in lattice simulations.^ If the chiral 
restoration turns out not to be a second order phase transition, but rather a cross-over 
or weak first order transition (finite but large correlation length), as the lattice data 
seem to indicate,a we expect the same qualitative features, though less pronounced. 

In summary, we used OPE and RGE to derive exact sum rules at finite (T, p) in 
asymptotically free theories. These exact sum rules strongly constrain the qualitative 
shape of SF's, especially near phase transition region. We urge whoever parameterizes 
a SF, e.g. in the QCD sum rule calculations or to interpret lattice simulations, to 
incorporate these exact constraints. In the future, we plan to generalize these results to 
baryonic currents and analyze their phenomenological consequences in greater detail. 

We thank Profs. Fried and Miiller for organizing a very enjoyable workshop. This 
work is supported in part by the US Department of Energy. 
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